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We present a quantum multi-modal treatment describing Electromagnetically Induced Trans- 
parency (EIT) as a mechanism for storing continuous variable quantum information in light fields. 
Taking into account the atomic noise and decoherences of realistic experiments, we model numer- 
ically the propagation, storage, and readout of signals contained in the sideband amplitude and 
phase quadratures of a light pulse. An analytical treatment of the effects predicted by this more 
sophisticated model is then presented. Finally, we use quantum information benchmarks to examine 
the properties of the EIT-based memory and show the parameters needed to operate beyond the 
quantum limit. 
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One of the steps towards the realization of quantum 
computation is a device that allows the coherent stor- 
age of information. The Heisenberg Uncertainty Princi- 
ple (HUP) sets a limit on the quality of stored informa- 
tion that depends on direct measurement and subsequent 
reconstruction. Much experimental and theoretical re- 
search is directed towards quantum memories for light 
to circumvent this classical benchmark. To realize such 
memories, methods that provide a coherent interface be- 
tween large atomic ensembles and light fields have been 
proposed. 

A scheme using the off-resonant interaction of a light 
field with a large ensemble of three level atoms was pre- 
sented in [1] . Off-resonant Faraday rotation was also used 
as a mechanism for mapping quantum states of light onto 
atoms [2 . The storage of a light field was shown to be 
possible by controlling the spatial distribution of atomic 
shifts in optically thick ensembles of three level [3", "T and 
two level atoms [5 . Probably the most actively studied 
technique to achieve a quantum memory for light utilizes 
Electromagnetically Induced Transparency (EIT) [B^/T . 

Experiments using EIT in atoms carried out in a 
sodium magneto-optical trap (MOT) [6 and in hot ru- 
bidium vapor cells [8 have demonstrated the storage of a 
light pulse for a few milliseconds. In solid state systems a 
storage of more than Is has been achieved using photon 
echo techniques [9 . The quantum nature of single pho- 
ton Fock states was shown to be preserved when stored 
and released from a MOT [lOl [TTJ [12], and theoretical 
studies have proposed methods to enhance the storage 
efficiency in these experiments [13l [Ml [15] . 

Although controlled storage of the amplitude and 
phase quadratures of a probe field at the quantum limit 
using EIT has not been achieved, experimental efforts in 
this direction have shown the delay of the two quadra- 
tures of a continuous wave beam [16^ and the transmis- 
sion and delay of vacuum squeezing [171 HI]- Improve- 



ments in the efficacy of this system can be made with a 
better understanding of the sources of excess noise and 
loss. The transfer of the sideband statistics from op- 
tical fields to atoms also requires further investigation. 
In this work we develop a model describing the storage 
of the signals contained in the sideband amplitude and 
phase quadratures of a light pulse in the presence of deco- 
herences and associated atomic noise, and use quantum 
information benchmarks to show the quantum nature of 
the transfer. 

In the first part we present theoretical models that de- 
scribe the multi-mode propagation of an amplitude and 
phase modulated pulse and the storage of its informa- 
tion onto atomic states in EIT-based memories. A nu- 
merical phase space treatment of light storage treats sev- 
eral sources of inefficiency present in current experiments. 
Linearized Maxwell-Bloch equations are then solved an- 
alytically in the weak probe approximation to explain 
the behaviour of the atomic noise and to give an expres- 
sion for the time-bandwidth product of this system in the 
presence of decoherence and finite atomic density. 

Next, we develop criteria that quantify parameters for 
which EIT based memories are able to store information 
in the quantum regime. Several criteria have been de- 
veloped in the past to distinguish classical and quantum 
distributions of states in other quantum information pro- 
tocols, such as teleportation or quantum cryptography. 
Signal-transfer coefficients T, and conditional variances 
Fcv, have been used as a state independent measure to 
analyze the effectiveness of teleportation experiments in 
the presence of non- unity gain [191 HOI [H] . We propose to 
implement the TV diagram to define benchmarks for the 
storage of continuous variable information and identify 
the parameters required to enable a transfer of informa- 
tion that outperforms any classical strategy. 



I. MODEL 
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Previous theoretical work has characterized the effi- 
ciency of EIT as a delay line for continuous variable quan- 
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turn states [T4^,'22^. Considering a three level atom, under 
conditions where there is a pure dephasing rate between 
the ground states, information can be slowed down within 
a narrow frequency window, and no additional noise is 
introduced beyond that which is necessary to preserve 
the canonical commutation relation of the field [22 • The 
width of the transparency window depends on the cou- 
pling beam power and the atomic density. Controlling 
the coupling beam in time allows storage of the informa- 
tion within the atomic sample. 

This storage process can be understood as follows. The 
coupling beam prepares the atoms initially in state |1) 
through optical pumping. When a weak probe propa- 
gates in the medium under EIT conditions, coherences 
are created between the two ground states of the atoms. 
These coherences arise from a quantum interference be- 
tween the two fields and acquire the sideband informa- 
tion of the probe pulse during its compression inside the 
medium. After the compression, most of the probe field 
energy has been transferred to the coupling beam and 
left the cell at the speed of light. At this point in time 
the atoms possess the frequency information of the probe 
within the transparency window, distributed in momen- 
tum space. When the coupling beam is turned off, the re- 
maining energy in the probe field leaves the medium with- 
out affecting the information stored. The information 
will be saved provided the readout is performed before 
the decoherence processes have affected the atomic state. 
When the coupling beam is switched back on, the probe 
beam is regenerated with the supply of photons from the 
coupling beam and leaves the medium while reading the 
spin state of the atoms. The main constraints are that 
the signal has to be encoded at frequencies within the 
transparency window and that the compressed pulse has 
to fit within the sample size. When these criteria are 
satisfied, the efficiency of this process is close to unity. 

Experimental investigation of this effect requires opti- 
cal sources at or below the shot noise limit, which is only 
possible at some modulation frequency around a carrier. 
Modeling this spatio-temporal quantum information ac- 
curately therefore requires a model which contains the 
quantum state of a large number of modes of the light, 
which we provide in this paper. We solve this problem 
numerically and then analytically to calculate the degra- 
dation of the signal and added noise during the stor- 
age process in the presence of decoherence mechanisms. 
Specifically, we consider dephasing affecting the ground 
state coherence and also allowing an exchange of popu- 
lation between the two ground states. 

We approximate the atomic structure by the three level 
atomic Lambda system shown in Fig. [l] where the two 
atomic ground states are degenerate and the transitions 
are addressed experimentally with orthogonal circular 
polarizations. The switching of the coupling beam can 
be done adiabatically or abruptly if the pulse is totally 
compressed within the medium [6] [23l |24] . Here we only 
consider the case where the coupling beam is switched 
abruptly, although other theoretical works show that the 



way the coupling beam is shaped in time enhances the 
efficiency [T5]. 

In the following calculation, we consider the simultane- 
ous storage of both quadratures of the probe when am- 
plitude and phase modulations are encoded within the 
EIT bandwidth. The preparation of the state can be 
achieved experimentally by passing a light pulse through 
amplitude and phase modulators sequentially. Provided 
the modulation frequency is larger than the Fourier width 
Auj of the pulse, classical information is encoded onto its 
sideband oo at the shot noise limit. 

We are interested in the envelope of the probe field, 
so the problem will be solved using the rotating wave 
approximation. The envelope operator is denoted by 
and its commutator is [£{z^t)^ £{z\t^)] = ^S{t — 
zj c—(t' —z' I c)\ where L is the quantization length, taken 
to be the length of the cell and c is the speed of light. Let 
X^(ijj) be the quadrature operators of the input probe 
field at the sideband frequency uj. For the amplitude and 
phase quadratures we have X^^J^uj) = fin(^) + ^/n(~^) 
and X~^{ij) = -i{Sin{^) - S-^{-uj)), respectively. 

We can write X^{lu) = 2a^{ij) + SX^{u), where 
(^^{uj) is the coherent amplitude encoded onto the probe 
via optical modulation, and 6X^(00) its quantum fluc- 
tuations. The power spectral density S^{uj) of a signal 
is the Fourier transform of the autocorrelation function 
and obeys the relation |22l |25] 

S^{co)5{co+co') = Ux^{uj)X^{uj')). (1) 

When normalized to the detection bandwidth, chosen to 
be much smaller than the applied modulation frequency, 
the measured power spectrum is 

5±H = |{|1±H|2). (2) 

We also need to define the noise floor of the signal by 
introducing the measured power spectrum without signal 

y±H = |(|,5x^HH. (3) 

For the input probe state we then have S^{u) = 
4f (a^^(cj))^ + V-^{u;). The signal is defined as 
4|;(a^^(cj))^ and the noise as ]^^(cc;), which is unity for 
a shot noise limited laser beam. 

If this state is inefliciently stored with some frequency 
and quadrature dependant linear loss r]^{uj) and if some 
excess noise with amplitude ^oise(^) generated by the 
memory, we have 

5±t = r,±(c.)5± + 1 - rj^iu;) + K± .J^)- (4) 

The term 1 — r]^{uj) corresponds to uncorrelated vac- 
uum noise, common to any system in the presence of 
linear loss. This noise term preserves the purity of a 
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quantum state and is necessary to preserve the commu- 
tation relations of the output state. On the other hand, 
a device generating some excess noise V^^^^{lj), trans- 
forms an initial coherent state, V^'^{uj) = 1 into a mixed 
state where V^^{u) = 1 + ^noise(^) > 1- In the following 
section we will calculate r]^{uj) and ^noise(^) using phase 
space simulations in the positive P representation. 



A. Stochastic Simulations 

We treat the probe beam as a general quantized field 
with longitudinal spatial dependence and the coupling 
beam flc{t) as a classical field. The atoms are all pre- 
pared in state |1) before the probe enters the cell via 
optical pumping induced by l^c- In this study we as- 
sume the coupling beam Rabi frequency to be 10^ times 
larger than that of the probe, ensuring that no atoms 
will move into state |2) due to optical pumping induced 
by the probe and the coupling beam will not be depleted 
throughout the storage process. Its dependence on z will 
therefore be ignored in this treatment. The validity of 
this approximation in the presence of decoherence is dis- 
cussed in the next section. When both beams are of 
comparable strength and are both treated as quantum 
fields, a strong correlation also takes place between them 
[26l [27] . This correlation will not affect the quantum 
statistics of the probe in our study so the coupling beam 
can be treated as a classical field. 




S{z,t) 



7o,7. 



FIG. 1: EIT level structure. S{z,t) is the envelope operator of 
the probe field, and ^c{t) the coupling beam Rabi frequency. 
Almost all the atoms are pumped into state |1) initially. 7 is 
the spontaneous emission rate from the upper state and 70 , jc 
are mean decoherence rates between the two ground states for 
pure dephasing and population exchange respectively. These 
two quantities are usually referred to as I/T2 and 1/Ti in the 
field of magnetic resonance. 



The master equation of this system is 



(5) 



where p is the reduced density matrix of the field and 
atomic variables and Hint is the interaction Hamilto- 
nian. We define locally averaged atomic dipole operators 
cFij{z^t) for the \i) — \ j) transition given by [7| |22] 



aij{z,t) 



nA5z ^ 

Zk^OZ 
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where A is the cross-sectional area of the beam, n the 
atomic density and 5z an infinitesimal slice of the medium 
containing N atoms. In the rotating wave approxima- 
tion, the interaction Hamiltonian of the Lambda system 
is then 



NTi 

— [gcF^i{z, t)S{z, t)^Qc{t)d-s2{z, t)^H.c.]dz 



(7) 

where g is the coupling strength on the probe transi- 
tion. The Cis are Liouvillians modeling the decays due 
to spontaneous emission from the upper state |3), and 
are defined by 



Zk^Sz 



1 



P^t^^li) (8) 



where for simplicity we assume the decay rates 7 from 
the upper state to be the same for both transitions. 

^[1^2]^ accounts for an off-diagonal dephasing rate 70 af- 
fecting the ground state coherence and arises from elastic 
collisions or atoms moving in and out of the interaction 
region defined by the probe beam quantized mode. Its 
expression is 



^[1,2] yP\ = 70 2^ (c^iiPCTn - -(JiiCJiiP - ^P^ii^ii) 

Zk^5z 

+70 (^22/5*^22 - \^l2^t2p-\p^l2^22)- (9) 

Zk^5z 

This term only describes the situation where the 
atomic population remains in state |1) during the whole 
process. If the pumping preparation is not optimum or 
if inelastic collisions are non- negligible, a population ex- 
change term C^i2] ^^^^s to be introduced. 

It is defined as 

rcoll \ ;(\ ^, \ ^ { ±.k c:±.k ^ :i.k ±.k c: ^ c:±.k ±.k \ 
^[1,2][P\ =^c (C^12P^21 - 2^12^2lP - 2^^12^2l) 

Zk^Sz 

+7c ^^2lP^l2 - ^<^21'^12P- ^P<^21<7i2)(10) 

Zk^5z 

^[12] [p\ ^^^^ affects the off-diagonal terms in the den- 
sity matrix in the same way as C^i'^^[p\^ but as the 
sources of these two decoherence processes are different 
we monitor them separately. It should be noted that this 
last term does not account for a pure loss of atoms out 
of the system, due to possible atomic motion out of the 
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interaction region or atoms moving into other hyperfine 
states. We also assume the mean dephasing rates de- 
scribing quantum jumps from |1) to |2) to be the same 
as the mean rates describing quantum jumps from |2) to 
|1) for simphcity. 



(a) 
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The normahy ordered characteristic function [29^ is then 




(b) 
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FIG. 2: Phase space numerical simulations of quantum infor- 
mation storage using EIT. Amplitude and phase modulations 
at 190 kHz are applied to the pulse. The decoherence rates are 
7o = 250 Hz, 7c = 100 Hz. a) 3D graph showing the storage 
of the probe amplitude quadrature on a time-space grid, b) 
are the variances of the input / output fields for the amplitude 
and phase quadratures, 1 corresponds to the quantum noise 
limit. These simulations used 2000 trajectories. 



To model this system, we used stochastic phase space 
methods, and worked with the Positive-P representation 
[28 . This phase space representation is computationally 
intensive but has the advantage of being exact as opposed 
to the truncated Wigner representation. We choose the 
following normal ordering of the operators 

(f ^(j|3,aj3,(j|2,<3-33,a22,<3-ii,(Ji3,a23,<3-i2,f ) (11) 

and define 



s(A,z)=n' 



.XiOiiz) 



(12) 



X(A,^)=Tr(p E{\,z)). 



(13) 



The equations of motion for xlA, ^) are calculated us- 
ing the master equation ([5| and the commutation prop- 
erties of the atomic and field operators. By taking the 
Fourier transform of the characteristic function equa- 
tions of motion, and assuming a large number of atoms 
in each slice Sz^ a Fokker-Plank equation can be found. 
The Stratonovitch corrections are small compared to all 
the other variables and are not included in the SDE. 
A set of nine c-numbers stochastic differential equations 
(SDE) describing the atomic dynamics with 18 uncorre- 
lated noise terms arising from atomic fiuctuations is then 
derived in the Ito form. Their expressions are given in 
the Appendix A. The Maxwell equations for the probe 
envelope in a moving frame at the speed of light, are 



-—a[z,t) Cr3(^,^) 

oz c 

--P[Z,t) (Tii{z,t) 

OZ c 



(14) 
(15) 



where Oi{z) refers to the ith operator in our normally 
ordered definition and A = (Ao...Ai...Aii) is a real vector. 



where the c-numbers a and P represent the operators 
and f , and (73, an correspond to the atomic operators 

^13, ^13- 

The evolution of a and /3 in space and time is computed 
when amplitude and phase modulations at a frequency 
0.0057 encoded onto a 5O/7 long coherent input state. 
The envelope of the field then presents two cycles in both 
quadratures and allows the information to be encoded at 
frequencies where the pulse is shot noise limited. We 
numerically evaluate the expectatio n values of the two 
quadrature operators (X+(z,t)) = a{z^t) -\- /3{z^t) and 
{X~{z^t)) = —i {a{z^t) — /3{z^t)) and their noise spec- 
trum S^{ij) = j^X^{z^ij)X'^{z^ —ij), where the averag- 
ing is done over a large number of trajectories in phase 
space. The noise floor V^{z^uj) is obtained by turning 
off the signal on the probe. We solved these stochastic 
equations using the numerical package XMDS [30 and 
chose parameters realistic to atom optics experiments 
with ^''Rb atoms. The atomic density was chosen to be 
10^^ cm^ with a total length of 12 cm. At the moment the 
pulse is inside the medium, the coupling beam is switched 
off abruptly and turned back on 5O/7 later. We chose for 
these particular simulations, a dephasing rate 70 = 250 
Hz and an inelastic scattering rate 7c = 100 Hz. 

Fig. [2] shows the results of this simulation where two 
quadratures of the multimode field have been stored in 
the presence of atomic noise. The stochastic simulations 
were averaged over 2000 trajectories. Fig. [2] (a) shows the 
propagation of the amplitude quadrature of the modu- 
lated pulse through the cell. The results are identical for 
the phase quadrature and are not shown here. We can 
see that the EIT-memory preserves the shape of the sig- 
nal with minimal distortion. To better quantify this we 
plot the power spectrum of the input and output fields 
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in Fig. [2] (b). The asymmetry in the transmission re- 
veals a frequency-dependent absorption of the pulse as it 
propagates through the system, characteristic of the EIT 
Lorentzian transmission window. We also see that 60 % 
of the classical signal is absorbed and that extra noise 
is added to the field. Using the previously defined nota- 
tion the transmission r]^{uj) = 0.40 and the excess noise 
^'oise(^) ~ 0.12. We will see in the last section if these 
conditions correspond to a quantum memory regime and 
describe the origin of the noise in the following section. 

It should be noted here that an iterative procedure was 
recently proposed to optimize the coupling beam shape 
and power [13^. Here, we chose our (time independent) 
coupling beam Rabi frequency by maximizing the output 
signal without decoherence, i.e we found a trade off be- 
tween off line center absorption and the compression of 
the pulse required to fit the sample. In this case the effi- 
ciency 77 was found to be 80 %, only limited by the lack of 
optical depth. Using the same procedure, at higher den- 
sities and re-optimizing the coupling beam strength, we 
found the transmission to be close to unity. Such time- 
bandwidth considerations are developed formally in next 
section B-2. 



B. Interpretation 

In this section we provide an explanation of the results 
found in the phase space simulations in the previous sec- 
tion. We first discuss the effects of decoherences on the 
losses and atomic noise introduced during the light prop- 
agation. We will show that excess noise can be under- 
stood as a preservation of the canonical commutation re- 
lations of the field in the presence of gain in the medium. 
We will quantify this by solving the Heisenberg-Langevin 
equations in the weak probe approximation in the case 
of information delay, and compare it with a more gen- 
eral theory of amplification and attenuation. We then 
describe the mapping and readout of the information en- 
coded on the probe, derive boundaries for optimum stor- 
age, and quantify the maximum information that can be 
stored in this system. As in our numerical simulations, 
the process will be solved when the coupling beam is 
switched off abruptly. 



i. The role of decoherences 

We will here focus on the noise properties of the EIT 
as a delay line to explain the excess noise observed. From 
the interaction Hamiltonian Eq. ([7|), we can obtain a set 
of Heisenberg-Langevin equations 

<3-i3 = -(7 + 7o/2 + 7c/2)(Ji3 + i^f((Jii -(J33) 
+ iQcO-12 + A3 



C^32 



C^12 
d 



dz 



S 



-(7 + 7o/2 + 7c/2)o-32 + i^ci^ss - ^22) 

i g'^£^cFi2 + F32 

-(70 + 7c)^i2 - W^cr^2 + i^lcFi^ + F12 
cri3 



(16) 



where we have included the decays of the atomic dipole 
operators, the sources of decoherence introduced previ- 
ously and their associated Langevin noise operators Aj 
describing the coupling of the atoms to vacuum modes of 
large reservoirs. The expressions for the Langevin cor- 
relations are calculated using the Einstein generalized 
equations 221 E] and the non-zero contributions are 
given in Appendix (B). The system of equations ([l6| will 
be solved to first order in f , 7c/7 and 70/7. To ensure 
efficient EIT, we will also assume \Vtc\^ ^ (770, 77c)- 

We first perform a steady state analysis of this system. 
Assuming the coupling beam Rabi frequency to be real, 
the atomic steady states are found to be 



(<3-ii) = 1 - 2— , (6-22) — , ((3-33) — 

7 7 7 

\ 9{^) \ ^9lQ,^\ \ ^7c 

{(^12) = ^ , (cri3) = -^(^) . (^23) = 

(17) 

We first note that the atoms are no longer fully 
pumped in the state |1) due to population exchange 7c 
and a non-zero dipole {0-23) therefore appears on the cou- 
pling beam transition. In this paper, however, we have 
assumed that the coupling beam is not depleted. In or- 
der for these solutions to be consistent, we then need to 
find the regimes where the coupling beam is negligibly 
absorbed. We do so by solving the following Maxwell 
equation for the coupling beam propagation 



dz c 



(^23), 



the solution for which is 



(18) 



(19) 



where d = ^ is the optical depth of the medium. Al- 
though the coupling beam intensity is absorbed linearly 
through the medium, a negligible depletion is guaranteed 
under the condition 



77c 



> 2d 



(20) 



which we will require in all the following calculations. 
We checked that this condition is verified in the above 
numerical analysis and the one presented in the last sec- 
tion. We also note that because of the pure dephasing 
70 a dipole (<Ji3) is created on the probe transition. A 
portion of the mean probe field is then absorbed by the 
medium by an amount e""^^, where ao = 
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We will now calculate the evolution of the probe quan- 
tum field as it propagates through the medium in the 
same approximate regime. To simplify the equations, 
the fast-decaying atomic variables will be adiabatically 
eliminated (^((Jis, 6-23) ^ 7), making these equations 
only valid over timescales larger than the spontaneous 
emission time, which is the regime of interest for EIT. 
We follow the same procedure as in [22 and solve the 
equations in the Fourier domain. Using the steady state 
solutions listed above, we can eliminate the second order 
terms in the probe power and negligible Langevin noise 
contributions using Appendix B. 

The Maxwell equation for the field amplitude quadra- 
tures can be solved to give 



gN . 

c 



-h{uj){z — s) 



-A{uj){z — s) 



^ - ijlo + 7c) 



(21) 



where fi{u;) = —io; (7 + 7(^/2); = 70 +7c is the total 
decoherence rate and the susceptibility of the medium is 
given by 



g'^N {jd - - (^33)) - i{^32)^c 

c m(^) 



(22) 



The first part of Equation ( [2T[ ) describes the absorption 
and phase shift of the probe propagating with a group ve- 
locity given by Vg = —uj/^{A{u;)) inside the EIT medium. 
The last two terms in Equation ( [21] ) correspond to atomic 
noise added to the field due to dephasing. 

We now need to calculate the power spectrum of the 
output state as a function as the input state using Equa- 
tions ( 1|21 ) and the Langevin correlations listed in Ap- 
pendix B. First we note that 

23?(A(cj)) 



(23) 



links the linear absorption with the atomic noise, a direct 
consequence of the fiuctuation dissipation theorem. This 
allows us to obtain 

S^{z,io) = 7j{z,Lo) S^Jlo) + (1 - v{z,L0)){l + Nf) (24) 

where 7]{z,lo) = e~2W(A(w))z^ ^^^^ 



Nf 



mFi2,Fl,])+comF^s,Fl,]) 
27on2+^2(2^ + ^p_3^^)- 



(25) 



We note that the noise power spectrum is phase inde- 
pendent, indicating that the response of the medium is 
the same for both quadratures of the field. As expected, 
the normally ordered Langevin correlations are respon- 
sible for excess noise on the output field, the strength 



of which depends on the ratio between normally ordered 
Langevin correlations and Langevin commutators. How- 
ever, because to first order in 7c/7, (F/3F13) = (see 
Appendix B), spontaneous emission is not responsible 
for excess noise. = 47c, so only the population 

shuffling term 7c transforms a initial coherent state into 
a mixed state. Using the notations defined previously 
we have ^noise(^) — ~ v{^))^f' The preservation of 
the commutation relations of the output field is ensured 
by the anti-normally ordered Langevin correlations that 
give the noise term 1 — rj. Indeed, When 7c = 0, we can 
check that we have have V^(z,cj) = 1. 

In order to understand why population exchange be- 
tween the ground states is responsible for noise, we will 
simplify the equations further and concentrate on side- 
band frequencies close to the carrier. 

We first solve for the steady states with the only source 
of decoherence being the shuffling terms 7c(^ii — <3"22) and 
7c (^22 — o'li). We find a new solution for the atomic po- 
larisation ((T13) and insert this into the Maxwell equation 



c 



This corre- 



to give {£{z)) = (fin)e^^ where a 
sponds to a population exchange driven amplification of 
the probe field inside the medium, the energy for which 
will be provided by the coupling beam, up to a limit also 
set by Equation (20). 



This shuffling term alone is, however, not physically re- 
alistic. As can be seen from the stochastic equation listed 
in the appendix and the Equations (16), the Liouvillian 



^[12] ^^^^ includes a ground state dephasing with mean 
rate 7c giving an extra linear loss <^c = ^ ^1 similar to 
ao. When solving for the steady state solutions in this 
case, with the same approximations as above, we find 
the net transmission close to zero frequency to be unity. 
The losses in fact exactly compensate for the gain, and 
the EIT medium no longer performs amplification. Even 
though the transmission that includes >C^{^2] ^^es not de- 
pend explicitly on 7c, this underlying gain term results 
in the creation of excess noise. 

Using Eq. (25) close to uo = and for 70 7^ we can 



find the noise to be 



Vr 



.7c 



noise — ^ (1 

70 



); 



(26) 



whereas for 70 = 0, Koise = 2az. 

A similar expression was found in the theory of two 
beam coupling developed in [34 . The presence of excess 
noise on the output field was interpreted from the theory 
phase insensitive quantum amplifiers. In [32 the signal 
to noise ratio of the optical field was shown to degrade in 
the presence of gain, and extra noise has to be inserted in 
the field equations to preserve the commutation relations. 
Precisely, it was shown that the output of an ideal linear 
amplifier with a gain factor G > 1, relates to the input 
state by 



(27) 



where Sl is a vacuum mode of the reservoir. The power 
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spectrum at the output of an ideal phase insensitive am- 
phfier is then given by S^^^ = GS^ + G — 1. 

We win now fohow the approach of Jeffers et al. [33] 
and win consider here a general theory for amplification 
and attenuation. By concatenating m amplifying and 
attenuating infinitesimal slices with linear amplification 
l-\-a6z and attenuation 1 — aSz, where dz = z/m we will 
calculate the noise properties of the field at the output 
of such a sample and compare it to the previous result 
based on the Heisenberg-Langevin equations. The power 
spectrum of the field at a slice m is given by 



9=* 



(1 



{a-a)z ^^^^^ 
m 



2a ^(1 

i=i 



r(5f„ -i) + i 
r 



(a-a)z J 



(28) 



By going to the limit m ^ oo, converting the discrete 
slices into a continuous array, we get when a a 



S^{z) = rj'iz) St + (1 - V'mi + N'f) (29) 



where 



TV} 



2a 



and r]'{z) = e 



{a—a)z 



(30) 



This general treatment allows us to assess the amount 
of excess noise present at the output of a system when 
gain and attenuation are known quantities. Eq. (26) can 



readily be found again by replacing a and a by their value 
in the EIT system close to zero frequency which validates 
this interpretation. 

We will now compare the signal to noise ratio and noise 
found by the present theory to the results given by the 
phase space treatment in the case of information delay. 
The signal to noise ratio for both quadratures is defined 

by 



4{a^{z)y 
VHz) 



(31) 



Fig. [3] shows the evolution of the noise {V{z) = 
1 + Koise(^)) and the signal to noise ratio as a func- 
tion of the depth of propagation in three different situ- 
ations. We consider the following decoherence combina- 
tions : (70, 7c) = (0,0.0057), corresponding to curve (i); 
ilOilc) = (0.0057,0.0057), the curve (ii) and (70,7c) = 
(0.0057,0), the curve (iii). Fig. [3]; a) and Fig. ||;b) are 
the noise results from the two approaches and Fig. [sj^c) 
and Fig. |3|d) are the signal to noise ratio results. 

Both theories are in good agreement. For curve (iii) 
there is no population exchange between the ground state 
and therefore the noise never exceeds the shot noise level. 
For curve (i), the noise increases linearly as predicted 
when 7o = 0. For curve (ii), the noise increases exponen- 
tially according to Equation (24). 

Fig. [sj^c) and Fig. [sjd) compare the signal to noise ra- 
tios results from the two approaches and again a good 



agreement is found between the numerical simulations 
and the analytical solutions. Even though the excess 
noise power is larger for (i) than for (ii) and (iii), the 
signal transmission is 100 % with 7c only, therefore the 
output signal to noise ratio is larger for (i') than for (ii') 
and (iii'). 



Phase space calculations 



Analytical calculations 



> 





FIG. 3: Signal to noise ratios and noise results for the numer- 
ical simulations : (a) and (c), and analytical solutions : (b) 
and (d). (i) and (i') corresponds to 7c = O.OO57 and 70 = 0; 
(ii) and (ii') corresponds to 7c = O.OO57 and 70 = O.OO57; (iii) 
and (iii') corresponds to 70 = O.OO57 and 7c = 0. The same 
parameters as for the phase space simulations of light storage 
were chosen. 



2. Light storage 

In this section we present an analytical model of the 
light storage protocol. Our treatment describes the 
transfer of information from the modulation sidebands 
of the probe beam to the collective atomic coherences in 
the sample and vice versa, taking into account the same 
decoherence effects and the finite optical depth. We will 
again consider fast switching and symmetric conditions 
for the writing and retrieval. 

For lossless information delay the probe pulse does not 
need to be completely inside the sample. However, light 
storage can only be performed efficiently if the group 
velocity is small enough such that the entire pulse is lo- 
cated inside the sample just before the switching of the 
coupling beam. In that case, information delay can be 
seen as light storage where the coupling beam has been 
switched off and back on immediately afterwards. Pro- 
vided no noise is introduced due to the switching, which 
we verified numerically, the excess noise on the output 
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field can be determined by the previous information de- 
lay study and the Langevin terms will be ignored in the 
following. 

The storage process is treated in three steps. First, 
we describe the mapping of sidebands of the pulse of 
duration T to the atomic coherences in momentum space, 
this is the writing stage. The second step, the storage 
time discusses the influence of the decoherences when the 
coupling beam is off. The last step, the reading stage^ 
is the mapping of the information stored in momentum 
space back to a probe field. We model the relaxation 
between the ground states with the decoherence terms 70 
and 7c introduced previously, in the same approximate 
regime. Similarly to [23 , to first order in f , two coupled 
linear equations can then be derived 



(^+r^) ^i2{z,t) = u£{z,t) 
where we introduce the quantities 

= 7d- 



(32) 
(33) 



d' = d 



(7 + 7d/2) 
7 + 7d/2 



(34) 



Vp describes the pumping rate of photons from the cou- 
pling beam and d' is the optical depth seen by the probe 
without coupling beam, in the presence of population 
shuffling. To simplify further the notations we also intro- 
duce 



X 



gN 



c 7 + 



o'ss) 



7 + 7d/2 



(35) 



Because these equations are linear, we will deal with 
the atomic and field as c-numbers and remove the hats 
on the operators. 

Writing stage 

We define the collective ground state coherence as the 
Fourier transform in space of the locally averaged ground 
state coherence operator cri2(^,t). 



cri2{k,t) 



LJo 



(36) 



We will see that under FIT conditions this quantity fluc- 
tuates with the same standard deviation as the input 
field. During the writing stage the state of the probe at 
each point in space can be found using Fqs. (32) and (33) 



in the frequency domain. As expected, the result is iden- 
tical to the deterministic part of Equation (21). We can 



then obtain the mapping of the field in co space to the co- 
herences in momentum space when integrating Fq. (33). 



We consider the memory to work in the linearly disper- 
sive regime, i.e the differential phase shift seen by all the 
spectral components of the field is the same. This allows 
us to change variables from coq to koVg when integrating 
Fq. (33) and to get 



ai2{k,t) 



dkoEin{{ko - k)vg)Vw{ko, t) (37) 



where Vw is a "distortion function" which quantifies the 
losses due to the finite bandwidth C{^) — 3?(A(a;)), and 
the finite length of the cell. It is given by 



Vw{kQ,t) 



. ^{iko-Ciik-ko)vg))L 



Au \ikoL — ({{k — ko)vg)L 

. ^{Tp-i{k-ko)vg)t _ 2 

^ V Tp-i(k- ko)vg 



(38) 



The integration of Fq. (37) is performed between k — 
Auj/{2vg) and k + Auj/{2vg) where Auj = 1/T. 

We now require the medium to be optically thick, d' 
/cL, and the frequency where the information is encoded 
to be smaller than the pumping rate Tp. Both these 
conditions ensure a high efficiency of the writing process 
as we will see. In this regime Fq. ([37[) reduces to 



ai2{k,t) 



^^9 
AuT, 



-(1 



koL, 



X / dkoSin{{ko - k)vg)smc{—-). (39) 



This equation describes a downsampling of the infor- 
mation from the probe field to the atoms due to a finite 
optical depth during its compression. The information 
is loaded at a rate Tp onto the collective ground state 
coherences. This process is much faster than the time it 
takes for the pulse to enter the sample (which is on the 
order of T) . When Fourier transforming back to the spa- 
tial coordinate at a time toff — we get an expression 
in the form of a convolution 



T^off 

^12 



{z) = — sinc(A^^) * [H{L)£,^{-z/vg)] (40) 
i p 



where H{L) is a top hat function defining the atomic 
sample boundaries. For the probe pulse to fit the atomic 
sample we then require the duration of the pulse to sat- 
isfy the relation T <^ L/vg. In this case, and provided 
the spatial components of the probe are well within the 
sample {Akz <C 1), there is no loss of information and 
Fq. (40) can be written 



-£in{-z/Vg). 



(41) 



In this process the statistics of the probe field is dis- 
tributed onto the atomic ground state as it propagates 
through the medium. Most of its energy is transferred 
to the coupling beam and has left the cell at the speed 
of light. 
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Storage time 

The coupling beam is now switched off at t = toff- The 
evolution of the atomic coherence and of the remaining 
probe inside the medium will again be solved in the adi- 
abatic limit. The system of equations describing this 
process is 



d_ 
dt 



= -(7 + 7rf/2)cri3 + (an - 
= -(7 + 7^/2)^732-^^*^^12 



(42) 



The Maxwell equation describing the time evolution of 
the spatial modes of the remaining probe inside the cell 
is then 



(^+i/cc+l/r)f(fc,t) =0 



(43) 



where 1/r = ^-^{{(Th) — (0-33)). After a time 3r the 
majority of the probe energy is then absorbed by the 
medium. The ground state coherences depend on the 
probe field to second order so this short process does not 
affect the stored information. 

Assuming that the coupling beam is switched back on 
at a time t = ton, we integrate Eq. (33) from t = toff to 
ton to obtain 



t 



'12 



-7d(^on 



f)^Wf 

^12 



ik) 



(44) 



which describes a simple exponential decay of the 
coherences over time due to a non-zero dephasing rate 

Reading stage 

To describe the reading stage we evaluate the coher- 
ences in the presence of a field on the probe transition 
in momentum space. We first solve for <Ji2(/c,t) inde- 
pendently of the probe field by combining Eq. (32) and 
Eq. (|33|. We obtain 



ar2ik,t) = e-^^'^'alT{k) 



where 



m = 



d/L — ik 



(45) 



(46) 



We then follow the same procedure as in the writing 
stage. The Maxwell equation ( 32 ) is solved in uj space to 
give 



I 



dko o-i°2"( 



- ^0 



)Vr{z,u;o) (47) 



where is a distortion function now affecting the trans- 
fer from the atomic coherences to the field and is given 

by 



'Dr{z,LOo) 



—(■ 



(ia;o-C'(- 



{iuj^T-C'{^))T) 



(d'/L- 



- 1> 



(48) 



where C{k) = ^{/3{k)). The integration of Eq. ( [47| ) is 
performed between uj — VgAk/2 and u + VgAk/2. Under 
the same condition as for the writing stage (high density 
and small enough pumping rate 7^) we obtain 



xL 



VgAkd' 



dkoa{n^^)sincA. (49) 



The downsampling also occurs when the information is 
transferred from the ground state coherences to the probe 
due to the finite optical depth. We can again transform 
this expression in time and space to obtain the field at 
the output of the sample 



£{z,t) 



d' 



-sinc(Acjt) * [H(T)cr^°2"(-Vgt)] (50) 



This expression can be simplified in the case where the 
compressed probe pulse fits entirely within the atomic 
sample, i.e when the duration of the pulse T satisfies the 
relation T <C L/vg. There is then no loss of information 
and Eq. (50) can be written 



S{z,t) 



(51) 



Using continuity arguments one can combine Eq. ([5l|, 
Eq. (41) and Eq. (44) one can then see that 



^XL 
Tr^d^' 



-'yditon—t, 



•off) C. 



(52) 



This expression relates the input and output states in 
the presence of pure dephasing and population exchange 
between the ground states. One can show that when 
= 0, the output is the perfect replica of the input 
state. 



The result given by Eq. ( 52 ) was obtained after making 
two main assumptions, L ^ VgT and Aco ^ Fp, which 
can be gathered in this inequality relation 



Vg/L < Acj < 



(53) 



The information has to be encoded at frequencies that 
satisfy this relation for perfect storage efficiency. The 
lower bound has to be satisfied for the pulse to fit the 
atomic sample. With a long input pulse, i.e a small spec- 
tral extent Acj, a high density or a weak coupling beam 
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is required, whereas for a short input pulse Vg/L can 
be made larger. On the other hand, the upper bound 
defines the minimum EIT bandwidth tolerable to mini- 
mize the losses. A short input pulse will require a large 
coupling beam power, whereas a weaker coupling beam 
power (narrower EIT bandwidth) can be used. The time- 
bandwidth product of the system (F^ xvg/L) can in fact 
be found to be d! ^ i.e the number of independent samples 
from the probe that can be faithfully stored depends only 
on the density. Not only does 7c introduce excess noise 
on the output probe mode, but it also reduces the time 
bandwidth product at a given density. At infinite density 
one can than store an infinitely broad probe spectrum. 

This result is similar to the one found in the general 
case of time varying coupling beam in [13 . With a cou- 
pling power calculated by equating the length of the pulse 
in the cell with the length of the cell, Q?^ = dj/T^ the con- 
dition Tdj 1 could be found. Our condition however, 
requires the full width at half maximum of the pulse to 
be smaller than the cell length. As also mentioned in [13] , 
the input and output pulses durations (which are identi- 
cal in our case) also have to satisfy the relation T^d <^ 1 
for the information to be imprinted onto the atoms be- 
fore the pulse is absorbed. We note that provided Eq. 
(53) is verified, the condition Q'^/dj <C is a sufficient 



condition for jdT ^ 1 to hold. This is the case in all the 
numerical simulations presented in this paper. 



II. QUANTUM INFORMATION BENCHMARKS 



a) 

|^)epr 
b) 



AM PM 

■e — B-^ 
■e — 



Memory 



Joint 
Measurement 



Classical Memory 








AM PM 




FIG. 4: (a) General schematics for characterizing an optical 
memory. A pair of EPR entangled beams are encoded with 
amplitude and phase quadrature information. One of these 
beams is injected into, stored and readout from the optical 
memory whilst the other is being propagated in free space. 
A joint measurement with appropriate delay is then used to 
measure the quantum correlations between the quadratures 
of the two beams, (b) A classical teleporter scheme used 
as an optical memory. The input state is measured jointly 
on both quadratures using two homodyne detection schemes. 
Analogous to classical teleportation the measured information 
is stored for time r before fed-forward onto an independent 
laser beam with a feedforward gain, g. The feedforward gain 
is analogous to a transmission of ^/^^{oJ) for EIT based mem- 
ories, (c) Quantum memory using EIT. The input state is 
stored in the long lived ground state coherence of three level 
atoms in A configuration. 



In this section, we will investigate the storage of opti- 
cal information from a quantum informatic perspective. 
We will benchmark the results obtained in our modeling 
against known quantum information criteria. These cri- 
teria, which are fidelity, signal transfer coefficients and 
conditional variances, will enable the determination of 
whether a quantum strategy has been used in the stor- 
age and readout of a quantum state; whether an EIT 
based quantum memory is possible in an experimentally 
realistic situation, and whether the output of the storage 
process is indeed the best clone of its input. 

Fig. [4] shows the schematics of our quantum memory 
benchmark. It was shown in Ref. [35] that the optimal 
classical measure and prepare strategy for optical mem- 
ory is the classical teleporter scheme as shown in Fig. |4] 
(b), and we therefore benchmark the performance of our 
EIT quantum memory against this setup. In this classi- 
cal scheme, the storage time can be arbitrarily long with- 
out additional degradation. However, two conjugate ob- 
servables cannot be simultaneously measured and stored 
without paying a quantum of duty [191 136]. Moreover, the 
encoding of information onto an independent beam using 
amplitude and phase modulators will also introduce an- 
other quantum of noise. In total, the entire process will 
incur an additional two quanta of noise onto the output 
optical state. 

Possibly the best known benchmark in quantum infor- 



mation protocols is the fidelity which measures the wave 
function overlap between the output and input states. It 
is given by 



^=(^in|Pout|^in). (54) 

which, in the Wigner representation can be written 
^ = 2^ // Win(X+,X-)Wout(X+,X-)dX+dX-. (55) 



For Gaussian states with coherent amplitude and 
power spectrum 5^, the Wigner functions is 



w(x+,x- 



2 (X+-2a+)^ (S~-2ce~)^ 



nS+S- 



(56) 



The fidelity of the classical teleporter scheme can be 
easily calculated [H] using formula ( 56|55 ) and gives 



2 e" 



V'(2 + Koise)(2 + K"oise) 



(57) 



where fc± = a± (1 -J?^)V(2 + Koise). Ktise are the noise 
variances of the output field for the amplitude and phase 
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quadratures and is the feedforward gain. For an ideal 
classical memory with unity gain, = 1^ and one can 
see from the previous argument that a coherent input 
state will give ^noise ~ ^' ilms giving a classical limit of 
^ > 0.5. It has been shown by Grosshans and Grangier 
[38] that when the fidelity of a teleporter JT > 2/3, the 
output state is guaranteed to be the best cloned copy of 
the input state. This fidelity limit called the no-cloning 
limit for teleportation corresponds to the addition of only 
one quantum of noise in the entire process. 



Cloning regime 



> 

■g 
L 



classical regime 




0.5 

Linear Losses 

FIG. 5: Fidelity as a function of memory loss for FIT for 
Oiin{oo) = 10,5,2,1. The non-classical and the no-cloning 
regimes are reached when > 1/2 and T > 2/3, respectively. 



The use of entanglement in the context of quantum 
teleportation, or for example EIT for quantum memories 
is necessary to break these limits. We now quantify EIT- 
based quantum memories using this criterion. There is a 
direct analogy between the feedforward gain, g and the 
EIT transmission y^r]{u). Fig. s] shows the behavior of 
^, as defined in Equation (57| with g = ^/r](uJ)^ with 
varying memory loss for different coherent state ampli- 
tudes. We note that the maximum amount of memory 
loss tolerable for beating both limits are dependent on 
the coherent amplitudes of the input states. This shows 
that fidelity is a state dependent measure. 

The formula for the fidelity can be extended to mixed 
input states using 



(58) 



Jeong et al. [40] showed that this formula can again 
be used to benchmark quantum information protocols. 
Nevertheless, characterizing quantum memory using the 
state dependent fidelity as a measure will be complicated 
for exotic mixed states. 

An alternative measure to fidelity for the characteri- 
zation of quantum information protocols was proposed 
by Grangier et al. [37] for quantum non-demolition mea- 
surement and by Ralph and Lam [19] for quantum tele- 
portation. This alternative uses the signal transfer coef- 
ficients, T^, and the input-output conditional variances, 
V^^ to establish the efficacy of a process. The conditional 
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FIG. 6: Diagram for total signal transfer coefficient T ver- 
sus conditional variance product V. The classical limit line 
shows the optimal performance of a classical teleporter. The 
passive loss limit defines the performance of an FIT-based 
memory that does not produce excess noise. The unity gain 
curve is obtained with increasing excess noise in an FIT sys- 
tem with no loss. The Amplification limit corresponds to the 
ideal amplifier limit. Regions A, B, C and D correspond to 
the quantum regime; Regions B and C represent the regime 
where FPR entanglement is preserved; Region D is the lossless 
amplification region; Region C denotes the no-cloning limit. 



variances and signal transfer coefficients are defined as 



T± = 



n. 



out 

± 

out 

±~ 

in 



n. 



(59) 
(60) 



where IZ^.,. 

out/m 



is the signal to noise ratio of the out- 
put/input field defined by 



'Mn/out 



^^^in/out-^ 



(61) 



/out 



We now define two parameters that take into account 
the performances of the system on both conjugate ob- 
servables 



(62) 
(63) 



Fig. |6] shows the plot of a TV-diagram. Similar to the 
fidelity, there are corresponding classical and no-cloning 
limits in the TV-diagram for a teleporter or an optical 
memory. It can be shown that a classical teleporter can- 
not overcome the T > 1 or F < 1 limits. By tuning 
the feedfoward gain, ^, a classical teleporter will perform 
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at best at the "classical limit" curve as shown in Fig. |6] 
Ref [19] shows that this classical limit can be surpassed 
using quantum resource (Region A). With limited quan- 
tum resource, it is possible to have an output state with 

V < 1 (Region B). When the input state is from a pair 
of entangled beams, this performance corresponds to the 
preservation of EPR entanglement at the output [39] . 
With a stronger quantum resource, T > 1 and V < 1 
can be satisfied simultaneously. Grosshans and Grangier 

showed that under these conditions the output state 
represents the best cloned copy of the input. The lower 
right quadrant of the TV-diagram (Region C) therefore 
corresponds to the no-cloning regime. 

We now characterize the EIT-based quantum memory 
in terms of the TV diagram. When an EIT system does 
not generate excess noise, the performance of the memory 
is described by the linear loss limit line. Assuming that 
the transmission through the EIT medium is identical for 
both quadratures, it can be shown that V = l — r]{u;) and 
T = 2ri{uo). We note that the result suggests an EIT with 
linear loss will surpass the classical limit independent of 
^Jr]{uo). This is because unlike the classical teleporter, 
the output state obtained from a linear loss EIT is not 
being measured throughout the transmission. Thus there 
is no measurement quantum duty for all transmittivities. 
Moreover, an input entangled state through a linear loss 
device will always preserve some entanglement at the out- 
put. 

However, when some excess noise is introduced in the 
storage process, T will decrease and V increase more 
rapidly. This is the case if for example some amplifi- 
cation is involved. Indeed we have seen in the previous 
section that G — 1 quanta of noise will be introduced 
for a lossless memory with a gain G in order to preserve 
the commutation relations at the output. The perfor- 
mance of the memory is then described by the ideal loss- 
less amplifier line on the TV diagram where one can show 
that V = G - 1 and T = 2G/(2G - 1). The optimum 
situation will then be when the gain of the amplifier is 
unity, so that T = 2 and F = 0. As the gain increases 
the memory no longer performs in the no-cloning regime 
and reaches region (D) where no quantum correlation ex- 
ists between the input and output states anymore, but 
the signal transfer is always larger than what a classical 
memory could achieve. 

There are other possible sources of noise that do not 
amplify the signal. For example, any transfer between 
the coupling beam and the probe via non-linear pro- 
cesses [23 or non-ideal polarizers will contribute to excess 
noise. We can introduce the excess noise phenomenolog- 
ically with V^oisei which we can assume to be quadra- 
ture independent. The TV performance is now given by 

V = l-7/(a;) + 14oise, andT = 2 7/(cj)/(l + 14oise). Unlike 
classical teleportation and in the ab sence of an amplifica- 
tion process, the EIT gain \/r]{u;) can only be less than 
or equ al to unity. If we assume perfect transmittivity 
with ^/r]{^J) = 1, increasing excess noise produces the 
unity gain curve in the TV diagram. We note that al- 



though input signal is perfectly transmitted, the excess 
noise leads to a degradation on both T and V. 



(a) 2 
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FIG. 7: Classical, EPR and no-cloning regimes plotted as a 
function of EIT linear loss and excess noise in (a) . A is the 
non-classical regime; B is the EPR regime and C is the no- 
cloning regime. These limits are drawn in fig (b) on a loss-gain 
plot and have been derived from Equation (29). Contrary to 
the case of (a), with a large enough gain and sufficiently low 
losses in the memory, Region (D) can be reached. 



We will now define the parameters required to reach 
the quantum regime in both cases. We plot these 
quantum regimes with excess noise versus linear loss in 
Fig. [rj^a), and gain versus loss in Fig. ^h). In particu- 
lar, to define Fig. ^h) we calculated the T and V cor- 
responding to the situation described by Equation (29). 



Then we found the gain and losses for which the EIT- 
performance crosses our benchmarks. These diagrams 
determine whether an experiment is sufficiently low noise 
and transmissive for quantum information storage. The 
no-cloning limit can only be surpassed when ^/r]{^J) > 0.5 
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FIG. 8: TV diagrams showing the performance of the EIT 
memory, (a) shows the evolution of its efficiency for three 
different 70 values, the dotted lines representing the loci for 
a constant 70 and varying 7c. (b) shows the evolution the 
EIT memory for three different 7c values, the dotted lines 
representing the loci for a constant 7c and varying 70. 



and Fnoise < 1 simultaneously in both cases. We note 
that similar figure of merit has recently been developed 
by Coudreau et al. [41 , during the course of this work. 

Using the numerical model presented in the earlier sec- 
tions, we investigate the parameters required to imple- 
ment the storage of optical information in the quantum 
regime. We model the situation where a fixed medium of 
length L = 12 cm with an atomic density of 10^^ cm^ is 
used to store an optical signal encoded on a pulse. The 
length of the pulse is chosen here to be 5O/7 and the in- 
formation is encoded on the quadrature amplitudes at a 
sideband frequency co = O.OO57. We plot the evolution 
of the EIT-performance as a function of the decoherence 
rates in Fig. |8] At zero decoherence rates, we note that 
(T, V) is not (2,0) because of the finite optical thickness of 
the EIT medium. With the above parameters the bounds 
in Equation ( [53| ) are satisfied by only one order of mag- 
nitude which makes the storage process non ideal even in 
the absence of dephasing. The no-cloning limit is how- 
ever still beaten in that case. We see that the evolution 
of the performances of the memory with 70 and 7c is 
radically different as predicted earlier. When the deco- 
herence rates increase, T monotonically decreases and at 
some decoherence rate value the quantum regime is not 



reached anymore in both cases. It is also important to 
mention that for any value of the couple (70, 7c) the re- 
gion D is never reached. The gain term from population 
exchange is always associated with loss so the EIT mem- 
ory can never get close to the lossless amplifier regime. 

We also wish to stress the difference between our figure 
of merit and the one used in [13 . The present figure of 
merit considers the signal to noise ratios for both quadra- 
tures at a given sideband frequency, as well as the amount 
of excess noise added to both quadratures of the field. 
The ratio between the total number of photons between 
the output and input also provides a figure of merit for a 
quantum memory that does not introduce uncorrelated 
extra photons in the light field output mode. 

III. CONCLUSION 

We have developed a quantum multi-mode treatment 
describing the storage of the quantum information en- 
coded on the sideband quadrature amplitudes of a light 
pulse using both stochastic simulations and an analyt- 
ical treatment for EIT. The two models included both 
the atomic noise and decoherence rates of realistic ex- 
periments. In our model we have, however, assumed 
an ideal three level atomic structure with incident light 
fields that have constant transverse spatial intensities, 
and a mono-kinetic atomic ensemble for which light is ex- 
actly tuned on resonance. We have also neglected the ef- 
fect of the back coupling from spontaneous emission into 
the light field, such as "radiation trapping" [43]. With 
these assumptions, the optimum sideband frequency for 
which the storage process can be efficiently performed 
depends mainly on the optical density and the coupling 
beam power chosen to set-up the EIT. We have also cal- 
culated the time-bandwidth product of the EIT mem- 
ory and have shown that it only depends on the atomic 
density. The dependance of the decoherence rates with 
density was not taken into account in this paper, but 
would therefore constitute an important figure of merit 
for quantum memories. 

We finally proposed the use of quantum information 
criteria to benchmark the performance of quantum mem- 
ories against an optimal classical measure-and-prepare 
scheme. We have shown that for typical decoherence 
rates in current experiments quantum information on 
the sideband quadrature amplitudes can be stored for 
milliseconds in the no-cloning regime, in the presence of 
small amounts of linear loss and excess noise. 
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V. APPENDIX A 

We list here the stochastic equations describing the 
evolution of the atomic c-numbers in the presence of 
a quantized probe field and a classical pump for one 
slice Sz of the medium. The c-numbers a and (3 repre- 
sent the operators E^{z^t) and E{z^t). The atomic vari- 
ables (cTs, (74, (75, (76, cry, (jg, (Jio, cTii) represent the oper- 
ators (6-13, 6-23, ^12, ^11, ^22, ^21, ^32, ^3i)- The cquatiou 



for (78 linearly depends on and (77 via the population 
preservation equation (76 + (77 + erg = 1 and is therefore 
not computed. The noise terms rij (for j = 1 to 18) are 
all delta-correlated and follow a Gaussian distribution, 
and have been normalized by -^^^ • The variables 70 , 

% and Ec are all normalized to the spontaneous emission 
rate 7. 



(74 



erg 



-(1 + 7o/2 + 7c/2)cr3 + Ec(T5 - a{l - 2gq - ^7) + y ^(^/^ ~ + ^(^/^ + ^3)^2 

-{acjA + Ec(J^){n^ - iua) + ;r^\/(7c + 7o/2)(l - c^r -cr6)(m7 + V^nio + mi2) 
+\/2(% + 7o/2)(74(ni4 - mis) 

-(1 + 7o/2 + 7c/2)(74 + a(79 + E^^CTq + 2(77 - 1) - 

+ (n3 + m4)/7 + (% + 7o/2)(l - cr6 - (77) (mis + m6)/4 



|-((74 - Ec/g)ni - ^((74 + Ec/g)n2 



-(7c + 7o)cr5 - Q^cTio - EcCTs 



-(75 (ni - m2) + {a{aQ - (77) + (7^ + jQ/2)as){ns - in^)/2 



^(^5 - ^^e) + 7r^\/Q^(7ii + (3(7^ + 1 - (76 - (77 + 27o(77 + 7c(^6 + CF7){^ri^ + V^^9 + ^^ll) 

2v27 2^ 



(76 = 1 - (76 - (77 - 7c((76 (77) a(7ii /5(73 



(ni + m2) - OLG^{n^ - in^)/2 



+ V'a(7ii + /5(73 + 1 - (76 - (77V'a(7ii + /5(73 + 1 - (76 - (77 + 7c((76 + Cr7)(n5 + m6) 
+ ^^(cr6 + 0-7)(n7 - ni2) - aCTii + ^0-3 + 1 - (76 - (77719 

+ V^aCTii + /5(73 + 1 - (76 - (77 V^Q^Cm + /^Cr3 + 1 - (76 - (77 + 7c((76 + (77) + 7o^7(^14 " mis) 



-/5(75(mi5 + ni6)/2 ^(ni8 - mi7) 

7 

(77 = 1 - (76 - (77 - %((77 (76) ^c(cr4 + CTiq) 



/2?7 



(ni + m2) + oi(JQ{n^ - in^)/2 



+a/2((^(7io + £^ccr3)(n5 + m6) + \/ ^ (^6 + cr7)(ni2 - u^) + -^^^ Ec{(Ja + ctiq) + 1 - (76 - (77(718 - nn) 



+\/2(/5(74 + ^ccrii)(ni4 - mi3) + /5(75(ni6 + mi5)/2 + Z^ (ni8 - mi7) 

V2^7 

-(7c + 7o)^9 - /^cr4 - ^cCTii + \J aaii + f3as + 1 - (76 - (77 + %((76 + (77) + 7o(77(-m8 + V^ng - inn) 

2V7 

+ ^;^^(^i4 + mis) + {P{(^6 - cr?) + (7c + 7o/2)c^ii)(mi5 + ni6)/2 - y^(7g(mi7 + ^is) 
<yio = -(1 + 7o/2 + 7c/2)^io + P(T5 + ^c(2cr7 + (76 - 1) + (7c + 7o/2)(l - cr6 - cr7)(n3 - m4)/4 + (ni6 - mi5)/7 
+^Y^(crio + Ec/g)nir + ^^(^lo - ^c/^)^i8 

^Tll = -(l+7o/2 + 7c/2)c^ll+Se(7g-/5(l- 2(76 -(77) + ^/2(7o/2 + 7J(7lo(n5+m^ 

+ ^r^\/(7c + 7o/2)(l - ere -cr7)(-m7 + V^nio mi2) 

7 
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27 



(64) 



VI. APPENDIX B 

In this Appendix, we list the non-zero Langevin corre- 
lations corresponding to the system of equations (16). 



(Fi3(^i, ^1)^13(^2,^^2)) 
x((7 + 7c + 7o)(<3-33) 
,^1)^^13(^2,^2)) 
X (27(6-33) 

(^13(^1.^1)^2(^2,^^2)) 



6{zi - Z2)S{ui + U2) 

nA 

27(^11) - 7c(<3-ii - 6-22)) 
5{zi - Z2)S{u;i + 002) 
nA 

2(7 + 70 + 7c) (^33)) 
d{zi — Z2)S{uJi + UJ2) 
nA 



(^12(^1.^1)^13(^2,^2)) = 



X (7c + 70) (^32) 

d{zi — ^2)^(^1 + ^2) 



(Fi2(^i, ^1)^12(^2, ^2)) 



nA 

X (7c + 70) (^23) 

_ S{zi - ^2)^(^1 + ^2) 

~ nA 



x((7 + 7c + 7o)(^33) + 7c(<3-22 +<3-ii) +27o(crii)) 

, . \\ _ ^(^1 - ^2)(^(^1 +^2) 



(^12(^1.^1)^12(^2,^^2)) 



nA 



X (7(^33) + 7c (<3-22 + <3-ii) +27o(cr22)). 



(65) 
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